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Abstract We investigate the organization of traffic flow on preexisting 
uni- and bidirectional ant trails. Our investigations comprise a theoretical 
as well as an empirical part. We propose minimal models of uni- and bi- 
directional traffic flow implemented as cellular automata. Using these mod- 
^ ' els, the spatio-temporal organization of ants on the trail is studied. Based 

on this, some unusual flow characteristics which differ from those known 
from other traffic systems, like vehicular traffic or pedestrians dynamics, 
are found. The theoretical investigations are supplemented by an empirical 
study of bidirectional traffic on a trail of Leptogenys processionalis. Finally, 
we discuss some plausible implications of our observations from the per- 
' spective of flow optimization. 

O ' 

> 

• rH 

X 

Traffic-like problems are relevant in various biological systems ranging from 
the motion of motor proteins in cells on a microscopic scale [14115] to the 
behavior of social insects on a macroscopic scale [4ll3l [TUIITT] . This has led to 
the development of various models that try to explain the collective behavior 
observed empirically in these systems. Besides their interesting transport 
properties, colonies of social insects can be seen as multi-agent systems, 
facing and solving various kinds of problems PfTOlITT] . For that reason the 
social insect metaphor has already been employed with great success in 
computer science [l][2] . This can also be extended to other fields, for example 
models of pedestrians dynamics where the concept of chemotaxis 15 , 13J has 
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been adopted for incorporating the mutual interactions of pedestrians [26 . 
A 'virtual chemotaxis' mechanism allows these models to reproduce the 
collective phenomena observed empirically in the dynamics of large crowds. 
In addition, also effects arising from the direct coupling of counterflowing 
streams of pedestrians or ants have been subject of extensive empirical as 
well as theoretical investigations p7 l [T0 i riT | fT9] . 

Taking into account these basic aspects, which will be addressed in more 
detail below, minimal particle-hopping models for traffic flow on preexisting 
uni- and bidirectional ant trails have been proposed [Tl lTTlfTB] . These models 
are implemented as cellular automata that have been previously used to 
successfully describe traffic flow in various other fields, e.g. physical systems 
[5] like surface growth and human transport systems on highways and in 
cities [251IB]. 

The main focus of the present work consists in the proposal of very basic 
and hence quite general models. We identify the main features of the flow- 
density relation of already established trails and present simulation results 
for the proposed models in the full parameter regime of attainable densities. 

In addition to the theoretical considerations, we present the first em- 
pirical traffic data from a bidirectional trail of Leptogenys processionalis. 
In order to facilitate the comparison with the predictions of the simpli- 
fied models, the proposed experimental setup, that is the choice of the ant 
species and of the type of trail, has to satisfy certain restrictions. In order to 
reduce complexity, we observed traffic flow of a monomorphic species on a 
preexisting trail. More details will be given in Sec. |4] and in the subsequent 
discussion. 

Nevertheless, the techniques employed might also apply in a more general 
setup. The data obtained this way are compared with previous empirical 
results {4ll3l[2T| and give some indications for more realistic extensions of 
the minimal models studied here which focus only on some key features 
of ant traffic (additional aspects related to non-equilibrium physics or the 
application to pedestrian dynamics are treated in [2UlfTOj ). 



2 The unidirectional model 

The model for unidirectional flow that we have introduced in [7] and that 
will be considered further in this paper can be seen as an extension of the 
so-called totally asymmetric simple exclusion process (TASEP) 28,[6], a 
well studied model in non-equilibrium physics. Ants will be modelled by 
particles that move stochastically along a one-dimensional lattice. In the 
following, we define the TASEP and then certain extensions are added that 
lead to the unidirectional ant trail model that we propose in this paper. In 
a further step the unidirectional model will be extended to the bidirectional 
one. 



1 In the following we will use the terms ant and particle interchangeably. 
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2.1 TASEP regime 

The TASEP is defined by a simple set of rules. Particles (representing the 
ants) hop in a fixed direction from one lattice-site i to the next one i + 
1. Each site can be occupied only by one particle (exclusion principle). 
Time evolution is continuous, which can be realized by a random-sequential 
update scheme: In each update step, one site i is chosen at random. If site i 
is occupied and i + 1 is empty, hopping takes places with probability q. If site 
i + 1 is blocked by another particle, nothing happens (exclusion principle) 
and another site is chosen randomly. A particle leaving the last site L will 
hop to site 1, i.e. motion takes places in a ring geometry. More realistic 
geometries are possible, but the main effects in our models are quite robust 
against a change of boundary conditions [22] . 

One important way of characterizing a traffic-like system is the funda- 
mental diagram [23]. It relates the average velocity V or the flow F with 
the density g of agents in the system. Both relations F(g) and V(g) are 
equivalent due to the hydrodynamic relation F — gV. Nevertheless, both 
will be used in the following since certain features can be seen more clearly 
using the average velocity, other using the flow. 

In the TASEP with random-sequential dynamics and N particles hop- 
ping on a lattice with L sites, flow F and average velocity V as a function 
of density g are exactly given by (see e.g. [28,6]) 

N 

V(g)=q(l-g), F(g)=qg(l-g) with g=-. (1) 

Mutual blocking of particles is the only mechanism of interaction, leading 
to a strictly monotonic decrease of average velocity with increasing den- 
sity. Flow and average velocity are directly linked to the spatio-temporal 
distribution of particles. The density g gives the probability of finding a 
particle at a site i, whereas 1 — g gives the probability of finding a par- 
ticular site being unoccupied. The distribution of particles is homogeneous 
on time-average. Due to increased mutual blocking the average velocity de- 
creases strictly monotonically with increasing density, which agrees with the 
behavior found for example in vehicular traffic [23 , 6 J . 

2.2 Cluster regime 

The model for unidirectional ant trails [7] is obtained by extending the 
TASEP. Now each lattice site can also be occupied by a pheromone-mark 
and/or one ant (see Fig.[T]left). If a pheromone mark but no ant is present at 
site i + an ant at site i will hop to site i + 1 with probability Q. If neither 
an ant nor a mark are present, hopping will take place with probability 
q < Q. Also pheromones evaporate: if one site is marked but no particle 
is present, this mark will be removed (evaporated) with probability /. If a 
particle is present and the site is chosen for updating, the mark will not 
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evaporate, reflecting the reinforcement of the pheromone-mark by the ants. 
The corresponding probability p(i, At) of finding a pheromone at site i is 
therefore given by 

, . . fl if i is occupied 

p(i,At) = l [{1 _ f)At (2) 

where At is the time-interval between the occupation of a site i by two suc- 
ceeding ants, called time-headway in traffic engineering. Obviously the time- 
headway determines the decay of the probability of finding a pheromone. 
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Fig. 1 Definition of the uni- (left) and bidirectional model (right): The sym- 
bols correspond to particles (representing ants) moving to the right (>), particles 
moving to the left (<), and pheromone marks (•). Although a natural trail is 
two-dimensional, the ant motion can be considered to be one-dimensional allow- 
ing to map the trail onto a one- dimensional lattice with L sites. Particles move 
stochastically along the lattice, with hopping rates q and Q. Free pheromones can 
evaporate with rate /. Java applets illustrating the dynamics of the models can 
be found at www.thp.uni-koeln.de/^aj. 

In the case / = (/ = 1) pheromones will never (immediately) evaporate 
and the TASEP-case with hopping probability Q (q) is recovered (Fig. |2J). 
More interesting properties arise in case of < / < I. The most surprising 
feature is the non-monotonic dependence of the average velocity on density 
for small evaporation rates / (Fig. [5] left). 

At low to intermediate densities, the average velocity stays constant 
(V = q). Beyond a certain threshold value a sharp increase of velocity can 
be seen. In the regime of high densities, the monotonic decrease known from 
the TASEP is found. Both regimes are a consequence of the incorporation of 
pheromone marks. Each ant is followed by a trace of marks (Fig. [3] middle). 
A succeeding ant perceiving this trace will hop with probability Q > q. If 
the preceding particle sees no pheromone mark, it will hop with probability 
q < Q. Then faster ant will catch up with the slower ones forming a moving 
cluster (see Fig. [3] middle and Fig. [3] right). Similar results are known from 
systems with particle-wise defects [TH] where each particle j has its own 
individual hopping probability qj. Here also clusters are formed and the 
average velocity stays constant v = q. An example from vehicular traffic 
is a platoon of cars following a slow truck. At very low evaporation rates 
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Fig. 2 Fundamental diagrams for the cellular automaton implementation of the 
unidirectional model: Q = 0.9, q = 0.2 and / = 0(o), 0.0002(a), 0.0008(V), 
0.002(«), 0.008(>), 0.02 (x), 0.08(o), 0.2(*), 1(«). In case of / = 1 and / = the 
average velocity V(g) and flow F(g) are exactly known from the TASEP. At low 
to intermediate densities, the average velocity stays constant for a suitable choice 
of evaporation probability /. 



or at high densities (small average distance or distance- headway) , all ants 
perceive pheromone marks. So on average their hopping probabilities all 
reach the same value Q, leading to TASEP-like features. The corresponding 
distribution of ants becomes homogeneous (see Fig. [3] left). 




Fig. 3 Space-time plots for the unidirectional model: Q — 0.9, q = 0.2, / = 0.001, 
g — 0.2. On the left the TASEP case (p = 0.9) is shown. Besides fluctuations, 
particles are distributed randomly. The other two plots show different stages of 
the cluster formation in the unidirectional ant trail model. Each ant (black) is 
followed by a trace of pheromones (grey). Due to evaporation this trace has a 
finite length. As indicated by the different slopes of the ants trajectories (middle) 
clusters move with different velocities depending on the presence of pheromones. 
Finally, in the stationary state only one moving cluster survives (right). 



2.3 A comparison with models of vehicular traffic 

As already mentioned, cellular automaton models are well established in 
case of vehicular traffic [5] ■ They seek to capture the complex behavior of a 
system composed of humans driving in their cars. Besides certain similari- 
ties, models of ant traffic should also reflect crucial differences to vehicular 
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traffic. One is the lack of some kind of velocity memory. In our model, ants 
reach their walking speed within one single update step. This is similar to 
the motion of pedestrians where the walking speed is also reached within a 
time less than 1 sec. In contrast, in vehicular traffic cars cannot accelerate 
instantaneously to the maximum velocity and so velocity is only increased 
gradually [SllESlE]. 

Furthermore, the behavior of ants can be expected to be much more 
homogeneous than that of drivers on a highway. Here a mixture of different 
vehicles with different maximum speeds, acceleration capabilities, and so 
on, but also of different kinds of drivers with different moods or attitudes, 
influences the overall behavior [23]. Therefore, the range of possible veloc- 
ities in vehicular traffic can be assumed to be larger than in ant traffic, 
at least in comparison with monomorphic species. As a consequence of the 
speed homogeneity, overtaking is rarely observed on ant trails at least under 
certain conditions [16] and has therefore not been included into our model. 

Also the definition of the update procedure itself can have a strong influ- 
ence on the dynamics. For modelling vehicular traffic a time-parallel update 
is widely used which performs a synchronous update of all lattice sites [6]. 
As only the occupation of the lattice at the moment of updating is used, a 
synchronous update incorporates time latencies of the drivers. But in case of 
ants the perception range is limited to their immediate environment [13,24. 
115] and the random-sequential update described above is more appropriate. 

The stochasticity in the model has mainly two sources. First, many in- 
fluencing factors are not known or difficult to include explicitly since the 
model would become too complicated. Thus they arc included in a statisti- 
cal sense through probabilities for a certain behavior. The second reason lies 
in the behavior of the ants themselves. They appear to possess an intrinsic 
stochasticity depending on the evolution of the particular species [13] . This 
is also widely used in applications [T|[2] . In vehicular traffic stochasticity is 
used to incorporate fluctuations in the driver's behavior which can lead to 
spontaneously formed phantom jams [25j . 

3 The bidirectional model 

For extending the unidirectional model to the case of bidirectional flow we 
have investigated several models [19,17 ■ As one common requirement they 
should reduce to the unidirectional model in case of vanishing counterflow. 
In the extension discussed here another lattice for ants moving in the oppo- 
site direction is added (see Fig.[T]right). Both lattices for ants obey the same 
set of dynamical rules. An extension to multiple lanes is then straightforward 
by adding further one-dimensional lattices. Experimental studies already 
addressed this issue j9][3] . Going further one could consider the extension to 
a two-dimensional model. Here the trail topology including the formation 
of lanes itself could be incorporated more naturally. Especially the relation 
between flow properties and trail topology could be investigated. Also the 
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analogies to already existing models of pedestrians dynamics surely provide 
a fruitful perspective [26|f27] . 

Overall, one deals with three hopping probabilities depending on the 
occupation of the nearest neighboring site in hopping direction. In absence 
of counterflow again the rates q and Q apply. In case of counterflow one 
additional rate K (see Fig. [T] right) is used. This is a crucial difference to 
most models of vehicular traffic which usually neglect the coupling of lanes 
in opposing directions 6J. Workers facing each other in opposite directions 
have to slow down due to information exchange or just in order to avoid 
collision 4,3,9,11]. Therefore, we choose hopping probabilities satisfying 
K <q < Q. 

A related problem is the organization of ant traffic at bottlenecks. For 
Lasius niger this has already been investigated experimentally |llj . In that 
case traffic flow is organized such that the number of encounters inside the 
bottleneck is reduced. This is quite different from our bidirectional model 
where the main effect originates from these head-on encounters (see Fig. [9] 
right). Here some similarities to models of pedestrians dynamics can be 
drawn [27,19]. Under crowded conditions separated lanes for each direction 
are formed dynamically. These lanes are stabilized by incorporating the 
desire of pedestrians to reduce the number of encounters with others in the 
counterflow. 

3.1 Equal densities 

In general, the fundamental diagram can be expected to depend on both 
densities pl and pn of left- and right-moving ants. That is, it has the form 
F = F(pl,Pr). For simplicity, we start the discussion with equal densities 
in both directions. Finally, we investigate the full fundamental diagram 
focusing on the generic properties and differences to the other cases. 

3.1.1 Cluster regime The average velocity roughly shows the behavior al- 
ready known from the unidirectional case (see Fig. [5] left), including the 
anomalous density dependence (sec Fig. [4] left). As the average lifetime of 
the pheromone marks is determined by the mean distance-headway (mean 
ant-ant distance) of ants in both directions, this regime only exists at very 
low densities. Similarly to the strictly unidirectional case, moving clusters 
for each direction are formed (see Fig. [5] left). With increasing density the 
lifetime of the pheromone marks increases such that they become present 
at every site and thus suppressing the mechanism of cluster formation. 

3.1.2 Plateau regime The generic property of this model is found in flow 
and cannot be observed directly in the average velocity [17,16 . For all 
values of /, flow roughly stays constant over a certain density regime (see 
Fig.[4]right). Due to the behavior of the average velocity one observes a shift 
of the beginning of the plateau in flow to lower densities with decreasing 
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Fig. 4 Fundamental diagrams for the cellular automaton implementation of the 
bidirectional model: Q = 0.9, q = 0.2, K = 0.1 and / = 0(o), 0.0002(A), 
0.0008(V), 0.002(<), 0.008(t>), 0.02 (x), 0.08(o), 0.2(*), 1(»). The average velocity 
shows the same non-monotonicity already observed in the unidirectional model. 
The main property of the bidirectional model is exhibited by flow. For intermedi- 
ate densities the flow is nearly independent of the density q and the evaporation 
rate /. 
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Fig. 5 Space-time plots for the cellular automaton implementation of the bidi- 
rectional model: Q = 0.9, q = 0.2, K = 0.1, / = 0.002. The left plot shows the 
formation of small moving clusters for Ql ~ Qr = 0.03. One observes one left- 
moving (grey) and two right-moving clusters (black) . With increasing density this 
regime vanishes as the pheromones are present at any site. Due to the coupling 
to counterflow a large localized cluster emerges (plotted for Ql = Qr = 0.2 in the 
middle and on the right). 



/. But at intermediate to high densities, there is obviously no dependence 
on /. This effect originates from the mutual hindrance by countcrflowing 
ants. Similar effects are known from systems with lattice-wise disorder |29j 
where the hopping probabilities depend on the actual position (e.g., due to 
construction works or car accidents in vehicular traffic) . In our case we find 
several localized clusters of ants in both directions (see Fig. fright). They 
form some kind of dynamically induced defects, as ants facing these clusters 
move with reduced hopping probability K. 



3.2 Different densities 

For natural trails, different densities pn and p L for each direction are more 
realistic. Ants moving back to the nest might be carrying load and there- 
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fore will behave quite differently from outbound ants pfl[9l[l2] . More complex 
models could incorporate this by also employing different hopping rates de- 
pending on the direction. But here we consider only symmetric hopping 
rates and focus on right-moving ants, treating the left-moving ones as coun- 
terflow. This is no restriction due to the symmetry 

F R (p L ,p R )=F L (p R ,p L ) (3) 




Fig. 6 Fundamental diagrams for the full bidirectional model. Average veloc- 
ity V R (pL, Pr) and flow F R (pL, Pr) are shown for the full range of densities 
(Pl,Pr) G [0, f] x [0, 1]. In case of vanishing counterflow {pL = 0) obviously the 
unidirectional case is recovered. 




Fig. 7 Extended fundamental diagrams for the full bidirectional model. Total 
flow F tot (pL, pr) and effective flow F cS (p^, Pr) are shown for the full range of 
densities (pl,Pr) € [0, f] x [0, 1]. Symmetries according to Q and ([5]) are found. 
For clarifying the generic features due to counterflow, the limiting cases pr = 
and pL = are not shown. 



3.2.1 Properties of the full fundamental diagram The behaviour of the av- 
erage velocity shows strong deviations from the characteristics of strictly 
unidirectional flow (see Fig. [5] left). Obviously, the regime of constant ve- 
locity only exists for a small area in the pl, p^-plane. Although counterflow 
leads to a reduced hopping rate K 7 the average velocity is even increased 
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by counterflow. As already observed for the special case of equal parti- 
cle numbers, a non-monotonicity is found at very low densities (pl,Pr) £ 
]0, 0.1] x]0, 0.1]. For larger densities, a strictly monotonic decrease is found. 
As a result, the average velocity in case of counterflow is lower than for 
strictly unidirectional flow (p^ = or p R = 0). Overall, the observed feature 
is rather an effect of the shared pheromone lattice than of the reduced hop- 
ping rate due to counterflow. As pheromones become present at nearly any 
site, no particle- wise disorder can be induced anymore. Therefore a quasi- 
TASEP case is attained in which the impact of counterflow can be neglected. 
With increasing density of counterflowing ants the impact of pheromones 
vanishes and hopping predominantly takes place with rate K . As already 
observed in the special case pl = p R , flow reaches a nearly constant value 
(see Fig. E] right). 

3.2.2 Additional quantities of interest Total and effective flow are non-zero 
quantities in case of different particle numbers or hopping rates for both 
directions. The total flow F tot (pL, p R ) measures the total number of moving 
ants per unit time and is independent of the direction of movement. For 
quantifying the effective number of ants moving in one particular direction 
the corresponding flow F cS (pl, Pr) is used: 

F tot ( PLlPR ) = F R ( PL ,p R ) + F L (p LlPR ) = F tot (p R ,p L ), (4) 
F cS ( P l,Pr) = F R (p L ,p R ) - F L ( PLlPR ) = -F^{ PR , PL ) . (5) 

As a consequence of employing the same set of hopping rates for both direc- 
tions, the quantities exhibit different kinds of symmetries due to ([3]). One 
also observes plateaus which are a consequence of the plateaus found in F R 
and F L for an intermediate range of densities (pl,Pr) (see Fig. [7|). Both 
features might be subject of first qualitative empirical investigations. The 
total flow might be of interest especially in the ecological context of col- 
lective transport like foraging [4,3 since it measures the total performance 
of the employed foraging strategy, in contrast to the effective flow. On the 
other hand, the latter allows to identify the overall direction of movement. 

4 A brief outlook to empiricism 

The discussed analogies between ant- and vehicular traffic can also be used 
for experimental investigations. Based on this we collected qualitative and 
quantitative traffic data |16| . As our models are quite simple, the experi- 
mental setup, as well as the observed species, have to meet some well defined 
requirements. The observed trail section has to exhibit a constant shape for 
the period of data collection which is the equivalent of a static road. To 
minimize complexity, intersections or branchings are excluded. These would 
be analogous to on- and off-ramps which are known to have a strong in- 
fluence on vehicular traffic [61123], As already pointed out (see Sec. 12.3ft the 
moving agents in both systems can be assumed to be quite different. For our 
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Fig. 8 Illustration of the experimental setup. The observed section A-B used 
for data extraction is located in the middle of the bidirectional trail. A narrow 
passage on the left leads to a huge lawn (framebox 1). On the right the entrance 
to the nesting site can be seen (framebox 2). 




Fig. 9 Illustration of the observed section. Ants following the established trail 
leave the observed section in the same order of entering (left). Therefore the n- 
th ant entering will also be the n-th one leaving. Counting is done separately for 
each direction. Ants facing each other in opposite directions slow down in so-called 
head-on encounters (right). 

investigations we chose ants of the monomorphic species Leptogenys proces- 
sionalis. The analogy in vehicular traffic would be a traffic flow composed 
out of many copies of the same car driven by nearly the same driver. This 
choice ensures that all moving agents are basically equivalent at least for 
each direction of movement. Especially no laden ants were observed. 



4.1 Methods 

The field studies were performed on the campus of the Indian Institute of 
Science, Bangalore, India. We obtained qualitative as well as quantitative 
data for one bidirectional single lane trail which was observed for about 
25 minutes. Similar investigations have been made for ten other trails of 
the same type which have been found to exhibit the same qualitative and 
quantitative properties. The criteria employed for the choice of the trails are 
for example the trail topology (single-lane, bidirectional) or the observed 
means of interaction (e.g., head-on encounters). Particular features will be 
addressed in more detail in the subsequent discussion. 
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4-. 1.1 Basic quantities and observations For collecting quantitative data 
we chose a relatively small section located in the middle of the trail (Fig. [8]) 
to reduce additional effects arising from the boundaries. For the same reason 
obstacles affecting flow are avoided. The trail itself connects the nesting site 
on the right to a lawn on the left. Regarding the ecological context, ants 
were not observed to carry any kind of load like nesting material or pieces 
of prey. Nevertheless the colony was still present for a few more days after 
data collection. 

Probably the most surprising feature is the apparent absence of overtak- 
ing. Although ants temporarily left the trail and were passed by succeeding 
ones, we never observed an ant speeding up in order to overtake. Making 
use of this observation one is able to extract traffic data based on so-called 
cumulative counting [23]. Therefore, we use video recordings of the observed 
section (Fig. [5]). Counting itself had to be carried out by hand as various 
video tracking systems failed. The time of entering and leaving of every sin- 
gle ant is measured. Each ant produces a datapoint for entering (t J + , n J + ) and 
leaving (tr_ , n J _ ) the observed section depending on the direction of move- 
ment j € {L, R}. The n-th ant enters the section at time t 3 + and leaves it at 
time t J _. n+(t) (n_(t)) denotes the number of ants that have entered (left) 
the measurement section up to time t. The curves defined by the datapoints 
(ij_, n J + ) and (tL, nL) are called arrival and departure functions respectively 



One directly observable quantity is the flow of ants which is given as the 
slope of the arrival or departure function (Fig. [ITJ right inset): 



Ati(n)=ti(n)-ti(n-l) and f ± ( n ) = —- (je{L,R}). (6) 



Here At±(n) denotes the time-headway introduced earlier. Under the as- 
sumption of approximately constant velocities v^in) (see ©) along the ob- 
served section the distance-headway can be calculated additionally: 



For determining density or average velocity, the length L of the observed 
section is needed. By tracing the path of one single ant and putting marks 
to a transparency on the video screen one obtains L in units of the body 
length of the ants. This method avoids errors arising from the observer's 
perspective or from biasing the trail by direct measurements. 

J^.1.2 Derived quantities As ants do not overtake, the n-th ant entering is 
also the n-th one leaving (Fig. [9] left). Therefore, the average velocity of the 
ra-th ant while passing the observed section of length L is given by (Fig. [10] 
left inset): 



(Fig. HOi. 





(AT R (n)) = t R {n) - t R {n) 



(v R (n)) 



L 



(8) 



{AT R {n)) ' 
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This method in particular makes use of the average velocity rather than of 
instantaneous measurements. So the measurement is less affected by fluctu- 
ations directly observable in the walking speed or direction of movement. 

To extract further information about the dependence of traffic flow on 
the spatial distribution of the ants it is necessary to measure density. Al- 
though we will not discuss the fundamental diagram, the method used here 
has already been employed successfully in case of strictly unidirectional 
traffic [16] . 

The instantaneous particle number for one particular direction j G 
{L, R} is given by 

N j (t) = n> + (t)-ri>_(t) = const, while t G [t{,t{ +1 [- 7 {j G {L,R}) . (9) 

Here t J k denotes the time of an entering or leaving event for one particu- 
lar direction j. In both cases the number of ants is changed by one unit. 
Based on this, one obtains the number of ants iV 3 ' (t) for each direction at a 
particular instant of time. Since this number is usually not constant during 
the time a specific ant needs to pass through the observation section, some 
averaging over the instantaneous particle number N J (t) becomes necessary. 

Due to the coupling to counterflow also the number of counterflowing 
ants has to be taken into account. Generally, the total time-evolution of the 
trail occupation is given by a list M of entering or leaving events for both 
directions: 

M = {U} with ti G {ffe,tfc} • (10) 

The elements of M are sorted by time. The time-average of the instanta- 
neous particle numbers affecting the n-th ant while passing the observed 
section (Fig. [TO] left inset) is given by 



fi<t*(ra) 



t i= t«(n) 
ti<t R (n) 

(N? f (n)) = {AT B R{n)) with B= N L (t t )(t l+1 -U). (12) 

Here the average particle number affecting a right moving ant is calculated 
and therefore the left moving ants are treated as counterflow (Fig. [TO] left 
inset). The instantaneous particle numbers for each direction are averaged 
from the time of entering t+(n) to the time of leaving t^(n) of the n-th 
right moving ant. To facilitate comparison with the proposed models, we 
introduce dimensionless densities by 

A") = ,„h, » and p? f (n) = L (13) 

In case of p c f(n) < 1 the travel time (z\T i? (n)) is assigned to the unidirec- 
tional case and to the bidirectional case otherwise. 
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Time t [sec] 



Fig. 10 Illustration ol data extraction. Entering and leaving n(t?) ants 

produce datapoints. Based on these, travel time (z\T fl (n)} and average velocity 
(v R (n)) are calculated. Also the instantaneous number of ants N R (t) within the 
observed section is extracted (see left inset). Based on time headway inflow /+ (n) 
is calculated (see right inset). Generally, datapoints for entering and leaving ants 
exhibit a pairwise structure (see middle inset). 

4-2 Empirical results 

We already made use of some of the qualitative results. Those are also 
reflected in the quantitative measurements validating the employed meth- 
ods. The first one is the apparent absence of overtaking. Therefore one 
finds a pairwise structure of datapoints of entering and leaving ants (see 
Fig. [10] middle inset). Overall, the departure function is roughly equivalent 
to the arrival function shifted by the average travel time. Obviously, the 
time headways are not changed much while passing the observed section. 
This indicates stable distance-headways (see[7]) which will also be subject of 
subsequent investigations. In case of strictly unidirectional traffic this fea- 
ture is a consequence of ants moving in platoons, and can easily be observed 
directly [Hj. 

4-2.1 Measuring flow The dominance of flow of right moving ants can be 
observed directly on the trail and is also indicated by the slope of the corre- 
sponding arrival and departure functions (see Fig. [JO]). For t 6 [0, 225] s flow 
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in both directions is nearly the same. Starting from t « 225 s the flow of 
right moving ants increases whereas the flow of left moving ants decreases. 
Starting from t « 1100 s, flows are again nearly the same. One behavioral 
pattern that affects data-collection directly are so-called U-turns. The av- 
erage rate for right moving ants is about 4% for the first two time-intervals 
and reaches 7% for the last one. Left moving ants performed U-turns con- 
stantly at a rate of about 7%. As flow in both directions is quite fluctuating 
this might indicate that U-turns are performed independently of flow itself. 
The existence of U-turns has also consequences for the data analysis since 
the pairwise structure of datapoints required by the introduced method of 
measuring average single-ant velocities has to be restored. 
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bidirectional 
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Fig. 11 Velocity distributions for a real bidirectional single lane trail. The veloc- 
ity of the right moving ants is hardly affected by counterflow. Average velocity as 
well as variance are nearly the same in both cases. Additionally, the mean value of 
left moving ants in absence of counterflow is nearly the same as for right moving 
ants. But in presence of counterflow a significant decrease is found. Generally one 
also observes a larger variance for left moving ants. 
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^.2.2 Velocity- and distance-headway distributions For the time-interval 
of data collection we extracted average velocities according to ((HJ) ■ In order 
to investigate the impact of counterflow we employ velocity distributions 
rather than fundamental diagrams. Although those can also be calculated 
using the described setup, a very large number of datapoints in comparison 
to the unidirectional case is needed. One would have to cover the whole 
Pl x /Oj?-pl ane instead of the one-dimensional p-line. As a first observation, 
one finds that the velocity distributions of the right moving ants exhibits a 
much lower variance than for the left moving ones (see Fig. [5] right). Also 
the average velocity shows no significant decrease in case of counterflow. 
This is quite different for the left moving ants. Obviously, the average ve- 
locity is significantly decreased by counterflow due to the right moving ants. 
Overall, one finds two velocities: the reduced one for left moving ants due 
to counterflow and the second one attained by all right moving ants and 
for the unidirectional case of the left moving ones. This can be explained 
by making use of the advantage of employing video recordings for data ex- 
traction. Head-on encounters are observed for both directions (see Fig. [9] 
right). But additionally the right moving ants take advantage of some kind 
of follow-the-leader behavior. They are occupying the center of the trail. As 
a consequence, left moving ants move off-center also in absence of counter- 
flow. Therefore the left moving ants are more affected by counterflow. But 
also generally variance is much larger for left moving ants. This might be 
caused by the lower density in this direction as a density dependence of the 
velocity distribution has been found for strictly unidirectional trails |16j . 
From the videos one observes that the movement of the left moving ants 
is slightly less directed than for the right moving ants. Therefore a more 
fluctuating movement is observed. This could be a consequence of the di- 
rect connection between trail pheromones and mass orientation known from 
many other species [TMT3] . 

Also the directly observable spatial pattern are affected as indicated 
by the distribution of distance-headways (see Fig. [T2|) . Distances between 
left moving ants are randomly distributed exhibiting a so-called negative 
exponential distribution. In traffic engineering ( 23J) this corresponds to a 
homogeneous spatial distribution of cars on the road. With respect to our 
models we therefore find the equivalent of the TASEP-regime (see Fig. [3] 
left). For the right moving ants generally an asymmetric distribution of 
distance-headways appears to be more appropriate. This indicates a devi- 
ation form the homogeneous spatial distribution found for the left moving 
ants. 



5 Summary and discussion 

We have introduced two minimal cellular automaton models for traffic on 
uni- and bi-directional ant trails. Similarities and differences between these 
models and the models for vehicular traffic or pedestrians dynamics have 
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right moving 

unidirectional bidirectional 




Fig. 12 Distance-headway distributions for a real bidirectional single lane trail. 
The left moving ants are distributed homogeneously along the trail. So distance- 
headways exhibit a negative-exponential distribution (grey). The right moving 
ants show a different spatial pattern. Obviously an asymmetric distribution of 
distance- headways (e.g. log-normal) is more appropriate (black). 



been emphasized. Especially for pedestrian dynamics there are several sim- 
ilarities which also are supported by recent experiments (see e.g. |27j and 
references therein) . The main interaction in both ant-traffic models is based 
on different kinds of dynamically induced disorder. In the unidirectional case 
disorder is assigned to the moving agents by a virtual pheromone field. Al- 
though this is quite speculative (as, to our knowledge, no empirical evidence 
for such an interaction has been found yet), the main feature of the model 
is quite robust. Probably most other mechanisms, which assign particle- 
wise disorder with respect to the hopping rates, would produce nearly the 
same pattern. For the bidirectional model dynamically induced effective 
site-wise disorder is the key mechanism. This mechanism is independent of 
the pheromone field as ants moving in counterdirection induce a change in 
the hopping rates (from q or Q to K). So this mechanism is based on the 
directly observable behavior exemplified in Sec. 14.21 

The stationary state in both models is characterized by the spatio- 
temporal distribution of ants on the trail. We have shown that the flow 
properties of the system are directly linked to empirically observable aggre- 
gation patterns. In particular, the bidirectional model shows a rich variety 
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of patterns. Based on the general properties of the bidirectional model, we 
have also pointed out qualitative features of additional quantities which 
might be of interest in case of bidirectional flow and could also be tested 
experimentally. 

In addition to the investigation of theoretical models, we have presented 
empirical data obtained by a simple experimental setup. First observations 
for one particular species Leptogenys processionalis at least qualitatively 
confirm the spatio-temporal patterns predicted by the models |16j (see 
Fig. [6|). In general, clustering phenomena seem to be a common feature 
in ant colonies [13] . On a quantitative level, velocity distributions and fun- 
damental diagrams have been obtained [4lll6j. Also here the main features 
of the unidirectional model have been found. For the full bidirectional case 
only few experimental data are available until now [2T|l4]. Although, in 
principle, the employed setup is capable of obtaining the full bidirectional 
fundamental diagram, more data are needed. Therefore, we have restricted 
the discussion to velocity- and distance-headway distributions for investi- 
gating the impact of counterflow. It turns out that in the case of strongly 
asymmetric flow, as observed here, ants move differently depending on the 
direction. With respect to our models we would need to introduce different 
hopping rates for each direction. Nevertheless, the mutual slowing down, as 
predicted by the models, could be confirmed. Also the impact of counterflow 
on the spatial distribution of ants along the trail was shown. 

The empirically as well as the theoretically observed organization of 
traffic flow raises fundamental questions on the advantage of such orga- 
nizations. Both cellular automata models produce moving clusters for an 
appropriate choice of parameters. In the unidirectional model this happens 
for a relatively wide density regime. The existence of these clusters corre- 
sponds to the regime of constant average velocity. In the TASEP the average 
velocity decreases strictly monotonically with increasing density as ants are 
distributed homogeneously along the trail. So moving in clusters enables the 
ants to keep on moving at a higher velocity in comparison with a homoge- 
neous distribution in TASEP. This is obviously achieved by reducing mutual 
blocking. One might argue that the fundamental diagram shows a maximum 
value of average velocity. This maximum is attained at the point of sharp 
increase from the cluster- to the homogeneous distribution. Nevertheless, 
a maximum of the average velocity at this point would be quite unstable. 
Even small fluctuations in density would lead to large fluctuations in the 
average velocity. In comparison, the cluster regime is quite stable against 
fluctuations in density. Therefore, in natural systems clustering is expected 
to lead to a decrease of travel time which can be seen as a user optimum. 

In the bidirectional case, moving clusters only occur at very low den- 
sities. The underlying mechanism is still the same as in the unidirectional 
model. At intermediate to high densities the main feature, namely the for- 
mation of localized clusters leading to a constant value of flow, emerges. 
This feature solely depends on the presence of counterflow. In this situa- 
tion flow is the crucial quantity. In an appropriate ecological context, the 
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outbound flow for example determines the number of ants travelling to the 
food source, whereas inbound flow is related to the amount of food carried 
back to the nest. Very different values of inbound and outbound flow would 
lead to a too large or too small number of ants at the source. So a vanish- 
ing effective flow, which is roughly independent of density fluctuations, also 
ensures a constant number of workers at a particular destination. As this 
feature involves flow, it can be interpreted as some kind of system optimum. 
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